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Abstract 

The purpose of this study was to develop a biobetter version of recombinant human interferon-p la (rhIFN-p la) to improve 
its biophysical properties, such as aggregation, production and stability, and pharmacokinetic properties without 
jeopardizing its activity. To achieve this, we introduced additional glycosylation into rhIFN-p la via site-directed 
mutagenesis. Glycoengineering of rhIFN-p la resulted in a new molecular entity, termed R27T, which was defined as a 
rhIFN-p mutein with two N-glycosylation sites at 80 th (original site) and at an additional 25 th amino acid due to a mutation 
of Thr for Arg at position 27 th of rhIFN-p la. Glycoengineering had no effect on rhIFN-p ligand-receptor binding, as no loss 
of specific activity was observed. R27T showed improved stability and had a reduced propensity for aggregation and an 
increased half-life. Therefore, hyperglycosylated rhIFN-p could be a biobetter version of rhIFN-p la with a potential for use 
as a drug against multiple sclerosis. 
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Introduction 

Multiple sclerosis (MS) is a chronic neurodegenerative disease 
affecting the brain and spinal cord, leading to symptoms, including 
blurred vision, muscle weakness, trouble with mobility and 
balance, cognitive and memory problems, and sensory distur- 
bances [1-3]. It typically occurs between the ages of 20 to 50, is 
more common in women than in men, and has a variable course 
[3,4]. As gold standard MS therapeutics, recombinant human 
interferon- P (rhIFN-P) products are widely used as a first-line 
treatment and have had a good long-term safety record over the 
last few decades [5,6]. Although rhIFN-P is currently competitive 
with newer oral medicines that provide improved compliance and 
tolerance, it is still uncertain whether injectable treatments can be 
completely replaced by oral drugs, at least until oral drugs 
demonstrate good long-term safety records [6,7]. This is partic- 
ularly important because drug safety is one of the biggest issues in 
MS therapeutics, as MS is not a life-threatening disease but a life- 
long disease (since it occurs at an early age) [4,6]. Therefore, the 
production of biobetter versions of rhIFN-P would fulfill consid- 
erable unmet needs in the MS therapy, both with respect to 



considerations within the pharmaceutical industry, such as 
biophysical stability and low production costs, and in medicine, 
such as fewer side effects, longer dosing intervals and route of 
administration. 

Protein modification with polyethylene glycol (PEG) or oligo- 
saccharide moieties is an approach often used to improve bioactive 
proteins, especially in regard to physical and thermal stability, 
increased solubility, protection against enzymatic digestion, 
increased circulating half-life, and in some cases, decreased 
immunogenicity [8-10]. Notably, next generation rhIFN-P 
therapeutics has been created using Fc fusions or PEGylation of 
rhIFN-P at its N- or C-termrnal region, or at cysteine residues. 
PEGylated-Avonex is currently undergoing phase III clinical trials 
[1 1]. However, production problems with rhIFN-P la still remain 
a concern for the development of PEG-rhIFN-P la, because of 
rhIFN-P la aggregation, and the increased production cost 
associated with PEGylation. 

Unlike PEGylation, glycoengineering requires no additional 
manipulation processes after the construction of a relevant cell line 
because glycosylation is a natural protein modification within 
mammalian cells. In addition, it is well known that glycosylation is 
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important for the activity of rhlFN-fi [12,13]. There are two types 
of therapeutic rhIFN-P in clinical use: rhIFN-P la (Avonex and 
Rebif), which is produced in CHO cells and is singly glycosylated, 
and rhIFN-P lb (Betaseron), which is produced in Escherichia colt 
and is not glycosylated [14—16]. In many reports, rhIFN-P la 
shows a higher specific activity, lower immunogenicity and 
decreased propensity for aggregation than rhIFN-P lb [13,16]. 
Runkel et al. demonstrated that this remarkable difference is due to 
glycosylation, in which the glycan moiety confers protein stability, 
solubility and biological activity. Thus, glycoengineering may be a 
promising approach for improving the biophysical properties of 
glycoproteins, such as their structural stability and pharmacoki- 
netics [17]. 

Although rhIFN-P la is remarkably more stable and active than 
rhIFN-P lb, manufacturing of rhIFN-P la still continues to suffer 
from variable levels of expression and stability in mammalian cell 
lines, which often results in a low product yield [18,19]. The main 
challenge for developing rhIFN-P la is aggregation and poor 
stability, which can occur at very different stages in the 
development processes, and can result in low host viability, low 
productivity and the development of precipitates in solutions 
[19,20]. 

The approach oudined in the present study used site-specific 
hyperglycosylation via site-directed mutagenesis, which resulted in 
the development of a new molecular entity, termed R27T. Site- 
specific hyperglycosylation of R27T was confirmed by western blot 
analysis, isoelectric focusing, enzyme immunoassay, PNGase 
treatment and quantification of sialic acid. R27T displayed 
superior stability, solubility, productivity and pharmacokinetic 
properties without loss of specific activity or alterations in ligand- 
receptor binding. 

Materials and Methods 

Ethics Statement 

All experimental procedures and protocols for animal study had 
been approved by the Institutional Animal Care and Use 
Committee of the Seoul National University (protocol #SNU- 
200909-33). All efforts were made to minimize suffering of 
animals. 

Gene Construction, Expression and Purification of rhlFN- 
(3 1a Glycosylation Analogs 

Glycosylation analogs were constructed by performing site- 
direct mutagenesis via PCR on wild-type human IFN-p. The full- 
length gene was recovered and cloned into the pMSG expression 
vector (Patent#. US20040038394 Al, PanGen Biotech Inc., 
Gyeonggi-do, Korea). Stable transfection into CHO cells was 
performed using a dihydrofolate reductase selection system and 
selected methotrexate resistance clones were grown in methotrex- 
ate selective medium. Purification of rhIFN-P la and rhIFN-P 
mutant proteins, such as Dl ION, R27T and R27T-(NITV)2 was 
performed by PanGen Biotech Inc., in the same way. Firstly, for 
their purification, culture fluid containing proteins was applied to a 
column of blue Sepharose 6FF (GE Healthcare, Buckinghamshire, 
UK), which was then eluted with 35% propylene glycol-based 
phosphate buffer. After elution, the eluate was sequentially loaded 
onto CM Sepharose FF (GE Healthcare) and to C4 Reverse 
phase-high performance liquid chromatography (RP-HPLC) 
(Vydac, CA, USA). Purified proteins were concentrated and 
diafiltrated using a 10 kDa cut-off membrane in a Millipore 
Labscale Tangential Flow Filtration System (Millipore, MA, USA). 
Finally, gel filtration chromatography was performed on a 
Sephacryl 100HR column at 2.5 mL/min. An additional product, 



termed R27TAGlyc, during R27T purification was also obtained 
from the C4 RP-HPLC purification step, which was separated 
from R27T in this step. R27TAGlyc consisted of singly or double 
glycosylated rhIFN-P mutein at a ratio of approximately 7:3 (data 
not shown). As control material, Rebif was purchased from Merck 
KGaA (Hesse, Germany), respectively. 

Analysis of Expressed Proteins 

Purified rhIFN-P mutant samples were assessed by SDS-PAGE 
and western blotting using Anti-rhIFN-P antibody (R&D Systems, 
MN, USA) as a primary antibody. To confirm the presence of 
glycosylation, 1 mg/mL R27T in 20 mM sodium phosphate 
monobasic dihydrate, pH 7.5, was treated with PNGase F (Sigma- 
Aldrich, MO, USA) at 37°C. Samples were taken at different times 
over a period of 1 hr. Deglycosylation was monitored by SDS- 
PAGE. Isoelectric focusing was used to determine the isoelec- 
tropoint (pi) of the protein and was performed using a pH 3-10 
isoelectric focusing (IEF) gel (Invitrogen, CA, USA). R27T in 
media samples was quantified using a rhIFN-P ELISA kit (IBL, 
Hamburg, Germany) according to the manufacturer's instructions. 

Glycosylation Site Confirmation 

LC/ESI MS/MS was used for site-specific glycosylation 
analysis of R27T. R27T was reduced and then digested with 
trypsin/ Glu-C. Presence of carbohydrate specific fragment ions, 
such as ml z 204 and 366, in the product ion spectra, was analyzed 
from glycopeptide ions. For identification of the exact site of 
glycosylation, PNGase F was used to remove N-linked glycans 
from trypsin/ Glu-C glycopeptides, and for conversion of Asn to 
Asp. Deglycosylated peptides were sequenced using LC /MS/MS. 

Monosaccharide and Sialic Acid Composition Analysis 

A Dionex HPAEC (ThermoFisher Scientific, MA, USA) was 
used to analyze monosaccharide and sialic acid composition. 
AminoTrap Columns (ThermoFisher Scientific) and Carbopac 
PA 10 analytical columns (ThermoFisher Scientific) were used at a 
flow rate of 1 mL/min at 30°C. Acidic sugars were analyzed by 
exposing the samples in 100 mM NaOH and 1 M NaOAc. All 
monosaccharides, including neutral and amino sugars, were 
analyzed by exposing the samples to 18 mM NaOH for 25 min. 
The waveform used for pulsed amperometric detection (PAD) was 
the Dionex default program for carbohydrates. 

Dynamic Light Scattering (DLS) and Atteuated Total 
Reflectance Fourier Transform Infrared (ATR-FTIR) 
Spectroscopy 

Hydrodynamic size was investigated using a Zetasizer Nano 
ZS90 (Malvern Instruments, Baden- Wurttemberg, Germany). The 
temperature in the Zetasizer chamber was equilibrated to 10°C. 
Each sample was measured in a disposable sizing cuvette (Sarstedt, 
Germany). Hydrodynamic size and polydispersity index (PDI) 
were calculated from the auto-correlation function using Zetasizer 
software, version 6.32 (Malvern Instruments). ATR-FTIR spectra 
(4000-600 cm" 1 ) were collected at 4 cm resolution using a 
Nicolet 6700 spectrophotometer (ThermoFisher Scientific) with a 
golden gate accessory (diamond crystal). The ot-helix, P-sheet, P- 
turn, and random coil contents of the proteins were estimated 
from the amide I region of ATR-FTIR spectra. Peaks of the amide 
I region were first treated by Fourier self-deconvolution and then 
curve-fitted using the Gauss and Lorentz formula with OMNIC 
Peak Resolve software (ThermoFisher Scientific). The area 
corresponding to each secondary structure was calculated accord- 
ingly and expressed as a percentage of the sum of areas. To 
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measure kinetic thermostability, rhIFN-P la and R27T were 
incubated over 96 hr at 37°C and decay curves were generated for 
each protein. Detectable rhIFN-P was quantified by using the 
cytopathic effect (CPE) assay at each sampling time. 

In vitro Antiviral, Anti-proliferative and 
Immunomodulatory Activity 

Antiviral activities were measured to determine the capacity of 
rhIFN-P to protect A549 cells against the CPE of a lytic virus over 
a range of rhIFN-P concentrations. World Health Organization 
natural rhIFN (NIBSC code: 00/572) was utilized as a standard. 
To measure antiviral activity, A549 cells were seeded in 96-well 
plates and serial dilutions of the proteins were added. Plates were 
incubated for 22 hr and encephalomyocarditis virus (EMCV, 1000 
TCID 50/mL) was added. Following a further 22 hr incubation, 
cells were dyed with crystal violet, at room temperature for 1 hr, 
and then the dye was extracted with 2-methoxyethanol. Absor- 
bance at 570 nm was then measured. For measurement of anti- 
proliferation effects, Daudi cells were seeded and serial dilutions of 
the proteins were added, followed by incubation for 48 hr. Cell 
proliferation assays were performed using an EZ-Cytox cell 
viability assay kit following the manufacturer's protocol. Immu- 
nomodulatory effects of drugs were measured by analyzing the 



presentation of MHC Class I in A549 cells. A549 cells were seeded 
in 100 mm dishes and treated with a serial dilution of each 
protein, followed by incubation for 48 hr. After harvesting cells 
and adjusting each sample for equal cell density with FACS buffer, 
FACS was used to measure the expression of MHC class I 
peptides. 

Molecular Modeling of the R27T/IFNAR2 Complex 

Molecular models of glycosylated rhIFN-P la were built from 
the crystal structure of wild-type rhIFN-P la (PDB ID: 1AU1) 
[21]. Mutation of arginine to threonine at the 27 th residue, and N- 
linked glycosylation of 1AU1 were performed using UCSF 
Chimera [22,23] and GLYCAM [24-26]. Due to the high amino 
acid sequence identity (about 30%) between IFN-P la and IFN-oc 
2a, initial IFN-P la/IFNAR2 complex structures were generated 
using structural alignment with a model of IFN-a 2a/IFNAR2, 
whose structure was previously determined by NMR based 
docking methods (PDB ID: 2HYM) [27]. For modeling N- 
glycosylation, one of the major oligosaccharide structures 
(FA2G2S2, F: Core fucosylated, A2: biantennary with both 
GlcNAcs as bl-2 linked, G2: two galactose linked beta 1-4 to 
antenna, S2: two sialic acids linked to galactose) was chosen using 
a structure obtained from the hydrophilic interaction liquid 
chromatography (HILIC) profiles of R27T. The oligosaccharide 
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Figure 1. Construction and testing of rhIFN-p glycosylation analogs. (A) A schematic of the rhIFN-p protein. Boxes represented the locations 
of the five a-helices. Each vertical line represented a position with potential for an additional N-linked glycosylation site, as predicted by NetNGIyc. 
Introduced N-linked glycosylation consensus sequence sites were showed by diamonds(?). Purified samples were separated by (B) SDS-PAGE and (C) 
western blot analysis. (D) Analogs were subjected to N-glycanase digestion for the indicated times. (E) IEF analysis was performed over a pH range of 
3-10. 

doi:1 0.1 371 /journal.pone.0096967.g001 
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Figure 2. LC/ESI MS/MS of trypsin/Glu-C digests of R27T and Rebif. (A) Extraction ion chromatogram at m/z 204 and 366 for product ion 
spectra at 21 min and 46 min, respectively. Amino acid sequencing of deglycosylated C3 and CI 3 peptides by LC/ESI MS/MS. (B) Deglycosylated C3 
peptide MS spectrum and fragmentation information for R27T (Lot: 12104DS01 and 12103DS01) and Rebif. (C) MS spectrum of the deglycosylated 
C13 peptide and fragmentation information for R27T (Lot: 12104DS01 and 12103DS01) and Rebif. 
doi:1 0.1 371 /journal.pone.0096967.g002 



was built using Carbohydrate Builder (Woods Group. 2005-2013, 
GLYCAM Web. Complex Carbohydrate Research Center, 
University of Georgia, Athens, GA., http://www.glycam.com) 
and was attached to the N-glycosylation site of the wild-type IFN-(3 
la/IFNAR2 structure. The final structures were minimized using 
Amber force field (AMBER99SB) [28]. During minimization, 
Amber parameters were used for standard residues, and the 
Antechamber module was used to make parameters for non- 
standard residues. Steepest minimization was performed (using 
100 steps) to relieve unfavorable clashes, followed by 100 steps of 
conjugate gradient minimization. Steepest descent and conjugate 
gradient minimization size was 0.02 A. Molecular visualization 
was done using UCSF Chimera and PyMOL (The PyMOL 
Molecular Graphics System, Version 1.2r3pre, Schrodinger, CA, 
USA). 

In vivo Pharmacokinetic Analysis in Rats 

Sprague Dawley rats were purchased from Orient Bio, Inc. 
(Gyeonggi-do, Korea), and maintained on a 12 hr light/dark cycle 
in a temperature- and humidity-controlled animal research 
facility. All rats were 8 weeks old, weighing 280-300 g at the 
beginning of each experiment. Each experiment included 9 groups 
(n = 3 per group). The femoral artery and vein were cannulated 
with a polyethylene tube (PE-50; Clay Adams, NJ, USA) under 
anesthetization. rhIFN-(3 substances (R27T, R27TAGlyc and 
Rebif) were bolus-injected intravenously (via the femoral vein; 



IV), subcutaneously (at the abdomen; SC), or intramuscularly (at 
the leg; IM) at a dose of 1 MlU/kg to rats (total injection volume 
of approximately 0.3 mL). Approximately 120 u\L of blood was 
collected via the femoral artery at 0 (to serve as a control), 1, 5, 15, 
30, 60, 120, 240, 480, 720, and 1440 min after IV administration 
and at 0 (to serve as a control), 15, 30, 45, 60, 120, 180, 240, 360, 
480, 720, and 1440 min after SC and IM administration. 
Approximately 0.3 mL of heparinized 0.9% NaCl-injectable 
solution (20 IU/mL) was used to flush each cannula immediately 
after blood sampling. Blood samples were centrifuged immediate- 
ly, and a 50 uL aliquot of each plasma sample was stored in a — 
80°C freezer. All animals were humanely sacrificed at the end of 
experiments by euthanasia method: C0 2 inhalation in a non- 
precharged 4L chamber of a moderate fill rate. Standard non- 
compartmental analysis methods were used to calculate pharma- 
cokinetic parameters (WinNonlin; standard version 3.1; Pharsight, 
CA, USA) [29,30]. 

Statistical Analysis 

A p-value less than 0.05 was considered to be statistically 
significant. Statistical analysis was performed using Statistical 
Package for the Social Sciences (SPSS, IBM Corporation, NY, 
USA). Statistical significance was determined by analysis of 
variance (ANOVA) with Duncan's multiple range test posteriori. 
All results were expressed as mean ± standard deviation, except 



Table 1. Sugar composition analysis of R27T and Rebif (mol/mol protein). 







Fucose" 


N-acetyl-glucosamine b 


N-acetyl-galactosamine b 


Galactose 3 


Mannose" 


Sialic acid' 


R27T 


2.4 


12.3 


N.D. 


7.7 


7.1 


3.2 


Rebif 


1.1 


4.7 


N.D. 


3.0 


3.5 


1.2 



a 2 M TFA for neutral sugars at 100°C at 4 hr. 

b 6 N HC1 for amino sugars at 100 a C at 4 hr. 

c 0.1 N HCI for sialic acids(NANA+NGNA) at 80'C at 1 hr. 

doi:1 0.1 371 /journal.pone.0096967.t001 
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Figure 3. Size distribution of rhIFN-p 1a and R27T. The analysis of (A) rhIFN-p 1a and (B) R27T was performed at a scattering angle of 90° and 
distributed sizes, PDI, and Zeta averages were shown, (r.nm: radius in nanometers). (C) rhlFN-(3s were quantified by CPE at each time up to 72 hr. 
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Table 2. Secondary structure ratios derived from the ATR-FTIR spectra of rhIFN-p la and R27T in solution. 



Samples Amide I region in ATR FTIR spectra 

a-helix (%) p-sheet (%) p-Turn (%) Random Coil (%) 

rhIFN-pia 17.3 37.0 24.5 21.2 

R27T 29.8 31.3 25.6 13.3 



doi:1 0.1 371 /journal.pone.0096967.t002 

for time to reach a C max (T max ), which was expressed as the 
median (ranges). 

Results 

Construction of the rhIFN-p Glycosylation Analogs, R27T 

To best maintain the structural and functional properties of the 
protein, additional N- glycosylation sites were not created in the 
five helical regions, consisting of amino acid residues A(2— 22), 
B(51-71), C(80-107), D(l 12-136) and E(139-162). A number of 
sites including Arg27, Asp39, Gln72, Asp73, Ser74, Ser75, Aspl 10 
and Glul37 were then screened for the additional introduction of 
a consensus N-glycosylation sequence (Asn-X-Ser/Thr, where X is 
any amino acid except Pro). Each site was analyzed using 
NetNGlyc software on the Center for Biological Sequence Analysis 
website (Technical University of Denmark) (Table SI). Based on 
our analysis, we selected the most flexible overhand AB loop and 
the highest mobility CD loop to minimize the structural 
modification, remaining at Arg27, Asp39 and Aspl 10 
(Figure 1A). A substitution at Arg27 of rhIFN-fi la was chosen 
as the most promising site, as it had the highest substitution 
potential score in NetNGlyc, whereas Asp 110 was chosen as a 
negative control. The substitutions made were Arg27 and Aspl 10 
with Thr and Asn, respectively. In addition, a hyper-glycosylated 
rhlFN-fJ model was created by extension of the sequence 
"ANITVNITV, termed (NITV) 2 " into the C-terminus. Both 
R27T and R27T-(NITV)2 showed a clear increase in molecular 
weight (by 26 and 30 kDa, respectively) in western blots 
(Figure IB). No change in molecular weight was observed for 
the D110N analog, compared with native rhIFN-P, at 22 kDa. 
This was confirmed by SDS-PAGE of rhIFN-P after purification 
(Figure 1C). Treatment with PNGase, which catalyzes the release 
of N-linked oligosaccharides, resulted in the appearance of three 
bands of 26, 22 and 18 kDa in the R27T sample (Figure ID). The 
upper and middle bands were doubly or singly glycosylated 
respectively, whereas the lower band contained only the peptide. 
Thus the molecular weight increment was attributed to additional 
glycosylation. Negative charge was increased in R27T and R27T- 
(NITV) 2 , which had pi values of 5.3-6.4 and 4.6-6.4 respectively, 
compared to native rhIFN-P, with a pi of 7-8 (Figure IE). 
Notably, R27T, which had one additional glycan at the 25 th 
amino acid, was selected as the lead protein because the amino 
acid change for the additional glycosylation site was minimal. This 
could be advantageous, as it could result in lower immunogenicity 
than that for R27T-(NITV) 2 . 

Glycosylation Site Confirmation of R27T Analogs with 
Additional Glycosylation 

Specific glycosylation sites were conclusively confirmed by 
determination of the carbohydrate attachment site on the 
polypeptide backbone of the protein. Glycopeptides were identi- 
fied by performing quadrupole-time of flight (Q;TOF) MS on 
trypsin/ Glu-C digests of R27T and the control Rebif, in which the 



C3 and C13 peptide fragment had potential glycosylation sites at 
Asn25 and Asn80 in R27T, and Asn80 in Rebif (Figure SI A and 
SIB). Glycopeptide ions were extracted and confirmed by 
examining for the presence of oxonium ions, such as ml z. 204 
(HexNAc) and 366 (HexHexNAc), with retention times of 2 1 min 
and 46 min, respectively, on the extracted ion chromatogram 
(Figure 2A). The presence of oxonium ions was again determined 
in the product ion spectrum at 21 min and 46 min (Figure 2A). 
The identification of N-glycosylation sites by MS depends on the 
specific deamidation of asparagines to aspartic acid within the 
consensus sequence NX(S/T) upon cleavage of the glycan moiety 
by PNGase F. The deglycosylated peptide was sequenced by LC/ 
MS/MS (Figure 2B and 2C). Glycosylation sites of R27T were 
identified at Asn25 and Asn80 in R27T. 

Monosaccharide and Sialic Acid Composition Analysis 

We performed glycosylation analysis to compare the relative 
proportion of monosaccharide content, mono-, di-, tri- and tetra- 
antennary structures with those of the reference, Rebif. Individual 
sugar residues were identified and quantitated as moles of 
monosaccharide per moles of protein. The monosaccharide 
contents in 1 mol of R2 7T and Rebif were shown in Table 1 . 
The contents of each monosaccharide and sialic acid per moles of 
protein were increased with additional glycosylation, compared to 
Rebif. All carbohydrate moieties consisted of fucose, N-acetylglu- 
cosamine, galactose, mannose and sialic acid without N-acetylga- 
lactosamine. Therefore, they were N-linked complex-type sugar 
chains. 

Biophysical Analysis of Protein Stability by DSC, DLS and 
ATR-FTIR 

DSC thermograms of rhIFN-P la and R27T in 20 mM acetate 
buffer pH 4.2 were evaluated to obtain thermal unfolding events 
(Figure S2). Unfolding transition temperature (T m ), calorimetric 
enthalpy (AH) and van't Hoff enthalpy (AH V ) of rhIFN-P la were 
61.90°C, 39.88 kcal/mol and 106.3 kcal/mol, respectively. In 
addition, T m , AH and AH V of R27T were 59.07°C, 36.87 kcal/ 
mol and 103.6 kcal/mol, respectively. Even though R27T 
exhibited lower T m than rhIFN-P la while having almost similar 
AH and AH V it did not exhibit any visible particles after 
withdrawn from the DSC scan (data not provided). DLS was used 
to observe hydrodynamic size of rhIFN-P la and R27T and 
existence of aggregates in the aqueous environment as well. DLS 
measurements of the proteins were provided in Figure 3. rhIFN-P 
la gave two volume distribution peaks at around 2.700 nm and 
9.825 nm with 98.9% and 1.3% volume ratio, respectively. 
However, R27T had only single size distribution peak at around 
3.722 nm by volume distribution. 

Polydispersity index (PDI) is dimensionless and scaled from 0 to 
1 where values less than 0.05 are rarely seen other than with highly 
monodisperse standards. However, values higher than 0.7 suggest 
that the sample has a very polydispersed distribution and is not 
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suitable for the DLS measurement. PDI values of rhlFN-p 1 la and 
R27T were 0.607 and 0.306, respectively. Since rhIFN-p la had 
higher PDI value than R2 7T, the protein was more polydispersed 
and might not be stable in the aqueous environment compared to 
R27T. Therefore, the result may indicate that rhIFN-(3 la possess 
instability issue in the aqueous solution compared to R27T. In 
order to investigate secondary structural stability of the proteins, 



ATR FT-IR was selected to analyze the amide group I (1700— 
1600 cm ) in proteins. The amide group I region of ATR-FTIR 
spectra was separated into nine peaks with Fourier self- 
deconvoluted spectra of rhIFN-(3 la and R27T. Each ratio of 
composite area represented the corresponding percentage of each 
structure from peak 1 to 9; peak #1, 2, 6 and 7 (1692 cm -1 , 
1676 cm , 1636 cm 1 and 1623 cm , respectively; |3-sheet), 
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Figure 5. Model of the N-glycosylated R27T/IFNAR2 complex. (A) front view and (B) side view. The complex structure oligosaccharides, R27T, 
and IFNAR2 were shown as orange, magenta, and cyan color, respectively. New hydrogen bonds between the oligosaccharide and complex were 
presented as red bold line in C (with IFNAR2) and D (R27T). Also, 27T and 80N in R27T protein were colored by yellow (A,B) and hydrogen-bonded 
residues in R27T and IFNAR2 were colored by green (C,D). 
doi:1 0.1 371 /journal.pone.0096967.g005 



peak #3 (1665 cm ; reverse turn), peak #4 (1655 cm ; ex- 
helix), peak #5 (1646 cm ; random coil), peak ^8 (1615 cm ; 
side chain vibration) and peak #9 (1598 cm" ; P-turn) were 
resolved accordingly. After resolving the peaks, relative percentage 
of the contents was calculated. Table 2 showed the relative ratio of 
a-helix, P-sheet, P-turn and random coil of rhIFN-P la and R27T. 
R27T contained more ot-helix and less P-sheet than rhIFN-P la. It 
was already known from the crystal structure of rhIFN-P la that 
there was no P-sheet secondary structure. However, recent 
spectroscopy studies give that increasing intermolecular P-sheet 
structure is a common feature of protein aggregation although 
interferon consists of only a-helix [31-33]. Since P-sheet contents 
of rhIFN-P la were about 5.69% higher than R27T, it might 
suggests higher potential of protein aggregations [32]. 

To confirm the actual effect of stability issue, the aggregation 
kinetics was shown in figure 3C. CPE assay was used to detect 
bioactivity of R27T and rhIFN-P la in accelerated condition 37°C 
up to for 96 hr. The bioactivity with aggregation of R27T 
compared with rhIFN-P la had a less decrease within the first 
96 hr. The calculated half-life for rhIFN-P la was 96 hr in 
contrast to 503 hr for R27T. 

Maintenance of In vitro Activity 

R27T and Rebif was analyzed for in vitro anti-viral activity 
together with an R27TAGlyc. Values of 299x1 0 6 , 252 xlO 6 and 
288xlCT IU/mg were obtained from the average activity of the 
three lots of R27T, two lots of R27TAGlyc and Rebif, respectively 
(Figure 4A). No decrease in antiviral activity was observed with 



additional glycosylation. We also examined the relative potencies 
of R27T, R27TAGlyc, and Rebif in anti-proliferative and 
immunomodulation assays. Similar activities were observed in 
these assays. In the antiproliferative assay, IC50s of 14, 22 and 
29pg/mL were obtained for R27T, R27TAGlyc and Rebif 
respectively (Figure 4B). In a FACS assay, which measured IFN- 
inducible expression of MHC class I on the surface of A549 cells, 
50% responses were observed at approximately 265, 359 and 
238 pg/mL, respectively (Figure 4C). 

Molecular Modeling of the R27T/IFNAR2 Complex 

An R27T structure was generated in silica, based on the crystal 
structure of human IFN-P la. The R27T/IFNAR2 docking 
structure was shown in Figure 5A and 5B. N-glycosylation was 
predicted to readily occur at the 80 th amino acid, regardless of 
oligosaccharide structure or the distance (solvent accessible surface 
area 76. 7A by Naccess V2. 1.1.1) at which glycosylation could have 
no effect on the R27T/IFNAR2 complex due to far distance 
between oligosaccharide and receptor [34]. By contrast, glycosyl- 
ation at the 25 th amino acid was more difficult to model in silico 
than at the 80 th residue, as the 25 th amino acid side chain was 
oriented towards the inside of the protein. However, "wet" 
experiments, including SDS-PAGE, confirmed that it was possible 
to glycosylate the 25 th amino acid, in a similar manner to that of 
the 80 th amino acid (Figure IB and 1C). A structure for R27T with 
the 25 th residue glycosylated was obtained from the in silico data. 
This showed that glycosylation of R27T at the 25 th amino acid 
could stabilize the interaction between R27T and IFNAR2. In 
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particular, new hydrogen bonds were formed between Thr44 and 
Asp51 of IFNAR2 and the oligosaccharide moiety of R27T, 
thereby increasing R27T/IFNAR2 complex stability (Figure 5C). 



Similarly, hydrogen bonding between Arg35 of R27T and the 
oligosaccharide moiety contributed towards the overall stability of 
the helical conformations within R27T (Figure 5D). 
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Table 3. Pharmacokinetic parameters for IFN-p after IV, SC, and IM administration of R27T, R27TAGIyc, and Rebif at a dose of 1 
MlU/kg to rats (n = 3). 




Parameters 


R27T 


R27TAGIyc 


Rebif 


IV 


AUC (xio 4 lU-min/mL) 


54.8±7.4* 


37.2±2.1 


40.3 ±2.4 


Un (min) 


88.0±4.4 


63.3 ±11.5 


61.3±21.6 


CL (mL/min/kg) 


1.9±0.2* 


2.7±0.1 


2.5±0.2 


V ss (mL/kg) 


90.1 ±8.7 


1 1 0± 1 7.1 


107 ±36.2 


SC 


AUC (xio 4 IUvnin/mL) 


12.2 ±1.5* 


4.1 ±1.3 


5.2±0.5 


ti/2 (min) 


423 ±154 


247±121 


180±69.3 


C max (lU/mL) 


235 ±54.0 


161 ±=14.9 


246±42.0 


Tmax (min) 


30 (15-45) 


30 (30) 


45 (30-45) 


F (%) 


22.3 


11.1 


13.0 


IM 


AUC (xlO 4 lU-min/mL) 


16.3±2.0* 


6.9±1.3 


7.3±2.1 


ti/2 (min) 


426 ±44.7 


534±211 


150±11.8* 


C max (lU/mL) 


374±59.7 


172±35.1* 


314±98.4 


Tmax (min) 


45 (30-60) 


15 (15-120) 


45 (30-60) 


F (%) 


29.7 


18.8 


18.0 



*Significantly different to the other groups (p<0.05). 
doi:1 0.1 371 /joumal.pone.0096967.t003 



In vivo Pharmacokinetic Studies in Rats 

Plasma concentration-time profiles for the rhIFN-P variants 
were obtained in rats after IV, SC and IM administration of 
R27T, R27TAGlyc and Rebif at a dose of 1 MlU/kg (Figure 6). 
Relevant pharmacokinetic parameters were listed in Table 3. After 
IV administration of the three rhIFN-Ps formulations, bi- 
exponential declines in plasma IFN concentration were observed. 
The area under curve (AUC) values for R27T were significantly 
higher than those of R27TAGlyc and Rebif by 1.47 and 1.36 fold, 
respectively. Thus, the total body clearance (CL) values were 
significandy lower in R27T than in the other rhIFN-P formula- 
tions. Moreover, the terminal half-life (ti/2) values of R27T tended 
to be higher than those of the other rhIFN-P formulations by 1.39 
and 1 .44 fold, respectively. The apparent volume of distribution at 
steady state (V ss ) values were comparable among the three rhlFN- 
P formulations. After SC administration, the AUC and extent of 
absolute bioavailability (F) values of R27T were significantly 
higher than those of R27TAGlyc by 2.95 and 2.01 fold, 
respectively, and those of Rebif by 2.33 and 1.72 fold, respectively. 
Moreover, the t 1/2 values of R27T tended to be higher than those 
of R27TAGlyc and Rebif by 1.71 and 2.35 fold, respectively. The 
peak plasma concentration (C max ) and T max values were compa- 
rable among the three rhIFN-P formulations. After IM adminis- 
tration, the AUC and F values of R27T were significantly higher 
than those of R27TAGlyc by 2.36 and 1.58 fold, respectively, and 
those of Rebif by 2.24 and 1.65 fold, respectively. The ti/ 2 values 
of R27T and R27TAGlyc were significantly higher than those of 
Rebif by 2.84 and 3.56 fold, respectively. The C m , LX values of 
R27TAGlyc were significantly lower than those of R27T and 
Rebif, and the T max values were comparable among the three 
rhIFN-P formulations. The ti/2 values after SC and IM 
administration were much higher than those after IV administra- 
tion, indicating that each substance is subject to flip-flop kinetics. 



Discussion 

We have constructed a biobetter version of rhIFN-P la, termed 
R27T, as evidenced the results of structure-function studies of 
rhIFN-P la and the type I IFN-P receptor. R27T contained a 
substitution of Thr for Arg at position 27 th of rhIFN-P la, which 
resulted in additional glycosylation at the 25 th position. Glycoen- 
gineering was successfully performed, as demonstrated by the 
bioassay, MS spectrometry and in vitro data. Notably, the in vitro 
assays, including measures of antiviral, anti-proliferation and 
immunomodulation, showed that additional glycosylation at the 
25 th amino acid (as a result of the R27T mutation) had no effect on 
ligand-receptor binding. However, these results were different to 
the results of Runkel et al [35,36]. Using alanine-scanning 
mutagenesis, they showed that the 27 th amino acid, arginine, 
was a solvent exposed residue important for antiviral and reporter 
gene activity. Although alanine-scanning mutagenesis is a widely 
used technique to determine the functional role of each amino 
acid, mutagenesis of the Arg27, which has a polar residue group, 
to alanine, which has nonpolar, aliphatic residue group, could 
cause different effects than mutation into threonine, which has a 
polar residue group. In addition, glycosylation at the 25 th amino 
acid could cause conformational changes without influencing 
receptor binding. 

Glycosylation could also affect molecular stability. Stability 
issues, such as aggregation, have been one of the biggest challenges 
for the production of therapeutic proteins. Moreover, these issues 
occur at very different stages in the development processes, 
including drug production, purification, storage and delivery. To 
address these issues, the behavior of rhIFN-P la and R27T in 
solution was investigated using biophysical methods. Conforma- 
tional stability (T m value) of R27T by glycoengineering compared 
to rhIFN-P la was observed by DSC. T m of R27T was less than 
that of rhIFN-P la. Although glycosylation generally enhance the 
conformational stability, the position of glycosylation on protein is 
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known to be a key factor to modulate the stability [37]. Thus, 
additional glycosylation at Asn25 on R27T seemed to have less 
effect on conformational stability. However, DLS and ATR-FTIR 
data suggested that rhIFN-P la was more susceptible to 
biophysical stability, such as aggregation, than R27T. In 
particular, ATR-FTIR indicated that a significant amount of Ot- 
helical structure had been retained, whereas intermolecular 13- 
sheet content (which reflects aggregation) was decreased by the 
additional glycosylation (at the 25' amino acid). These results 
were consistent with the DLS data since rhIFN-(3 la was more 
polydispersed than R27T. Furthermore, these results were also 
confirmed by the aggregation kinetic profiles, which demonstrated 
the beneficial effect of additional glycosylation on stability during 
accelerated storage conditions. This effect could be caused by 
enhanced biophysical stability of R27T by glycoengineering. In 
addition, R27T productivity was approximately 3-6 times higher 
than that of rhIFN-(3 la under low temperature conditions of 
30°C-34°C for 6-14 days without additional manipulations like 
low temperature adapted cells, using cytopore, intermediate 
temperature shift and perfusion for production (Table 4). It was 
likely that this was a consequence of the lower aggregation 
resulting from the additional glycosylation, since it was well known 
that the low productivity of rhLFN-p' la was caused by molecular 
aggregation. 

To investigate the effect of glycosylation at the 25 th amino acid 
on stability, a computational simulation of R27T/IFNAR2 
docking was performed. Intrinsic stability was improved by 
glycosylation at the 25 th amino acid, which enabled hydrogen 
bonding between residue 35R on R27T and the oligosaccharide. 
Structural simulations of R27T/IFNAR2 suggested that the 
oligosaccharide could also interact with IFNAR2, resulting in 
the stabilization of ligand-receptor binding interactions. These 
results imply that glycosylation may be important for fine-tuning 
glycoprotein bioactivity although X-ray crystallography experi- 
ments may be required to obtain a more detailed structure [38] . 

Glycans play a very prominent role in determining therapeutic 
efficacy, including the in vivo half-life. We therefore compared the 
pharmacokinetic properties of R27T, R27TAGlyc and Rebif. 
After IV administration, the pharmacokinetic parameters (CL, V ss , 
and t 1/2 ) of Rebif (1 MlU/kg) were consistent with those obtained 
in a previous study on the pharmacokinetics of native rhIFN-P, 
administered at a dose of 21 MlU/kg [39]. The AUC and CL 
values of R27T were significantly higher and lower, respectively, 
than those of Rebif and R27TAGlyc. Although litde information 
on the elimination route of rhIFN-P la is currently available, the 
reduced CL is likely due to the fact that R27T is a hyper- 
glycosylated form of rhIFN-p. The reduced CL may be attributed 
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Figure SI Amino acid sequence of Rebif and R27T. 

Shadow on amino acid showed predicted glycopeptides and circles 
showed potential glycosylation sites of (A) Rebif and (B) R27T. 
(TIF) 

Figure S2 DSC thermograms of rhIFN-P la and R27T. 

Measured T m values of (A) rhIFN-P la and (B) R27T was 
displayed above the peaks. Graphs also showed curve fitting of the 
DSC thermograms for rhIFN-P la and R27T, which was 
performed using a non-two-state model. 
(TIF) 
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